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Gonzalez  RR,  Cheuvront  SN,  Montain  SJ,  Goodman  DA, 
Blanchard  LA,  Berglund  LG,  Sawka  MN.  Expanded  predic¬ 
tion  equations  of  human  sweat  loss  and  water  needs.  J  Appl 
Physiol  107:  379-388,  2009.  First  published  April  30,  2009; 
doi:  10.1 152/japplphysiol. 00089. 2009. — The  Institute  of  Medicine 
expressed  a  need  for  improved  sweating  rate  (msw)  prediction 
models  that  calculate  hourly  and  daily  water  needs  based  on 
metabolic  rate,  clothing,  and  environment.  More  than  25  years  ago, 
the  original  Shapiro  prediction  equation  (OSE)  was  formulated  as  msw 
(g-m_2-h_I)  =  27.9-£'req'CEmax)~0'455,  where  £req  is  required  evap¬ 
orative  heat  loss  and  £max  is  maximum  evaporative  power  of  the 
environment;  OSE  was  developed  for  a  limited  set  of  environments, 
exposures  times,  and  clothing  systems.  Recent  evidence  shows  that 
OSE  often  overpredicts  fluid  needs.  Our  study  developed  a  corrected 
OSE  and  a  new  msw  prediction  equation  by  using  independent  data 
sets  from  a  wide  range  of  environmental  conditions,  metabolic  rates 
(rest  to  £450  W/m2),  and  variable  exercise  durations.  Whole  body 
sweat  losses  were  carefully  measured  in  101  volunteers  (80  males  and 
21  females;  >500  observations)  by  using  a  variety  of  metabolic  rates 
over  a  range  of  environmental  conditions  (ambient  temperature,  15- 
46°C;  water  vapor  pressure,  0.27-4.45  kPa;  wind  speed,  0.4-2.5  m/s), 
clothing,  and  equipment  combinations  and  durations  (2-8  h).  Data  are 
expressed  as  grams  per  square  meter  per  hour  and  were  analyzed 
using  fuzzy  piecewise  regression.  OSE  overpredicted  sweating  rates 
(P  <  0.003)  compared  with  observed  msw.  Both  the  correction 
equation  (OSEc),  msw  =  147-exp  (0.0012-OSE),  and  a  new  piece- 
wise  (PW)  equation,  msw  =  147  +  1.527-iireq  —  0.87-Z?max  were 
derived,  compared  with  OSE,  and  then  cross-validated  against  inde¬ 
pendent  data  (21  males  and  9  females;  >200  observations).  OSEc  and 
PW  were  more  accurate  predictors  of  sweating  rate  (58  and  65%  more 
accurate,  P  <  0.01)  and  produced  minimal  error  (standard  error 
estimate  <  100  g-m_2-h_1)  for  conditions  both  within  and  outside 
the  original  OSE  domain  of  validity.  The  new  equations  provide  for 
more  accurate  sweat  predictions  over  a  broader  range  of  conditions 
with  applications  to  public  health,  military,  occupational,  and  sports 
medicine  settings. 

thermoregulation;  modeling;  fluid  balance;  hydration;  fluid  replace¬ 
ment 


daily  water  needs  can  be  determined  from  “minimal”  water 
losses  and  expected  increases  in  different  water  flux  avenues 
(11,  22).  Metabolic  water  production  and  respiratory  losses 
often  offset  each  other,  and  fecal  losses  are  usually  small  (11, 
22).  Urinary  losses  primarily  depend  on  hydration  status  and 
osmolar  load,  but  sweat  represents  the  largest  potential  avenue 
of  body  water  loss.  Knowledge  of  sweat  losses  is  therefore 
critical  for  calculating  water  needs  for  active  populations, 
particularly  when  exposed  to  heat  stress  (11,  22).  For  example. 
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military  potable  water  planning  relies  on  water  tables  generated 
from  existing  prediction  equations  (15,  25),  and  similar  nomo¬ 
grams  have  been  generated  for  public  health  purposes  (11). 
Although  the  sports  medicine  community  recognizes  the  vari¬ 
ability  of  human  sweat  losses  and  the  need  for  individualized 
fluid  replacement  guidance  (23),  realistic  sweating  estimates 
for  sporting  scenarios  enable  the  same  water  planning  capabil¬ 
ity  as  in  military  logistics  and  facilitate  accuracy  of  phenom¬ 
enological  modeling  in  sport  (14).  There  also  are  applications 
for  occupational  settings  (12)  and  disaster  relief  efforts. 

The  Institute  of  Medicine  (IOM)  recently  set  U.S.  Dietary 
Reference  Intake  (DRI)  standards  for  water  and  electrolytes 
(11).  The  IOM  used  the  original  Shapiro  equation  (OSE)  (28) 
to  estimate  sweating  rates  and  calculate  daily  water  needs,  as 
well  as  daily  sodium  losses,  over  a  broad  range  of  activity 
levels  and  environmental  conditions.  The  OSE  was  selected 
because  it  was  the  best  equation  available  at  the  time,  but, 
because  of  its  limitations,  the  IOM  panel  determined  a  need  for 
improvements  in  the  “development  of  capabilities  to  predict 
hourly  and  daily  water  requirements  based  on  metabolic  rate, 
climatic  conditions,  and  clothing”  (11).  The  OSE  (28)  is  a 
sequel  to  the  Givoni-Goldman  model  to  predict  core  tempera¬ 
ture  (6).  It  was  developed  from  laboratory  experiments  on  men 
for  energy  expenditures  ranging  from  —15  W  (rest)  up  to  475 
W  (moderate  metabolic  intensity)  over  a  range  of  environmen¬ 
tal  conditions  [20-54°C  and  10-94%  relative  humidity  (RH)] 
while  wearing  shorts  and  a  T-shirt  or  what  is  now  obsolete 
military  clothing  and  equipment  (27,  28).  The  original  derived 
equation  is  shown  below: 

OSE;  sweating  rate  (g-m^-lT1)  =  27.9-£req-(£max)_0'455  (  7) 

where  £j-eq  is  the  evaporation  required  to  maintain  heat  balance 
at  any  given  core  temperature  and  £max  is  the  maximal  evap¬ 
orative  capacity  of  the  environment.  Equation  1  has  been  used 
widely  to  predict  water  needs,  assuming  the  fluid  intake  (1/h) 
replaces  the  expected  water  lost  by  sweating  in  a  heat-accli- 
mated  person  [sweating  rate  X  body  surface  area  (BSA)  X 
10^2,  1/h], 

Although  Eleq  is  a  final  outcome  solved  by  the  solution  of 
heat  balance  and  is  a  determinant  of  those  parameters  involved 
in  thermoregulation  (i.e.,  skin  and  core  temperature,  skin  wet- 
tedness,  and  heat  production),  this  variable  is  not  wholly  equal 
to  thermoregulatory  sweating.  Equivalence  is  based  on  the 
efficiency  of  the  sweat  secreted  to  cool  the  skin.  For  100% 
effectiveness,  all  sweat  must  be  evaporated  at  the  skin,  and  any 
variance  in  efficiency  or  heat  storage  or  imbalances  or  incon¬ 
sistencies  in  heat  exchange  properties  can  affect  the  sweat  rate 
equation  predictability. 

The  OSE  often  overpredicts  sweating  rates  for  conditions 
both  within  and  outside  the  original  experimental  £leq  and  Emax 
domains  of  validity  (3).  Recent  tests  of  its  robustness  during 
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more  prolonged  exercise  bouts  (>2  h),  higher  exercise  inten¬ 
sities,  cooler  temperatures,  or  modern  protective  clothing  in¬ 
dicate  poor  agreement  between  predicted  and  measured  sweat 
losses  (3).  Cheuvront  et  al.  (3)  found  that  the  overpredictions 
from  use  of  the  OSE  were  also  likely  associated  with  compli¬ 
cations  attributed  directly  to  non-sweat  losses  of  body  mass 
(NSL).  In  addition,  variability  in  estimations  of  dry  heat  losses 
(radiant,  R  +  convective,  C)  and  £max  from  a  clothed  individ¬ 
ual  can  lead  to  incorrect  calculations  whenever  imprecise  heat 
and  vapor  transfer  coefficients  through  clothing,  determined  by 
use  of  static  manikins,  are  subsequently  applied  to  dynamic 
conditions  based  on  variable  wind  and  walking  conditions  (1, 
7,  8,  10).  It  is  imperative  that  improved  equations  be  developed 
and  validated  that  predict  sweating  rates  over  wider  thermal 
environments,  with  and  without  outdoor  solar  loads,  higher 
metabolic  rates,  longer  work  durations,  and  with  contemporary 
military  clothing  including  the  new  combat  uniform  and  other 
modem  protective  uniforms  and  equipment  (11,  12,  15,  23,  25). 

The  purposes  of  this  study  were  7)  to  compare  the  accuracy 
of  OSE  with  measured  sweating  rates  during  extended  exercise 
(>2  h)  with  higher  metabolic  rates,  contemporary  clothing 
systems  including  modern  body  armor,  and  a  broader  range  of 
environmental  conditions;  2)  to  develop  a  new  equation  that 
more  accurately  predicts  sweating  rates  compared  with  OSE; 
and,  if  possible,  3)  to  develop  a  correction  to  OSE,  designated 
as  OSE,c,  that  is  applicable  and  easily  migrated  into  various 
existing  rational  and  operational  thermal  prediction  models. 
The  first  hypothesis  was  that  OSE  overestimates  sweating  rates 
over  extended  periods  of  exercise  and,  therefore,  leads  to 
excessive  estimates  of  water  needs.  The  second  hypothesis  was 
that  new  and  corrected  equations,  based  on  a  more  compre¬ 
hensive  database,  would  provide  a  more  accurate  estimation  of 
water  needs  to  encompass  most  military  (15,  25),  occupational 
(12),  and  public  health  scenarios  (11),  as  well  as  many  situa¬ 
tions  within  the  sports  medicine  community  (23). 

METHODS 

The  database  consisted  of  101  volunteer  subjects  (80  men  and  21 
women)  with  >500  observations  included  in  the  data  set  used  to 
develop  the  algorithms.  Experiments  were  conducted  for  the  purpose 
of  precisely  measuring  sweating  rates  over  a  broad  range  of  condi¬ 
tions.  Each  protocol  was  approved  by  the  appropriate  institutional 
human  use  review  boards,  and  all  volunteers  were  informed  both 
verbally  and  in  writing  of  the  objectives  and  procedures  of  the 
respective  study.  No  identifications  of  a  given  volunteer’s  personal 
records  were  present  in  the  spreadsheet  database. 

Raw  data  were  obtained  from  four  separate  environmental  chamber 
studies  and  one  field  study  conducted  at  the  U.S.  Army  Research 
Institute  of  Environmental  Medicine  (USARIEM)  and  from  one 
environmental  chamber  study  previously  conducted  at  Defence  R&D 
Canada,  Toronto  (DRDC)  (8).  Actual  sweat  losses  were  carefully 
measured  in  all  studies  by  weighing  volunteers  and  accounting  and 
correcting  for  non-sweat  losses  of  mass,  food,  and  fluid  intakes  (NSL) 
(2,  3).  In  brief,  sweat  losses  were  measured  by  accounting  for  the 
change  in  individual  nude  body  mass  (kg)  measured  on  a  electronic 
precision  balance  scale  (Toledo  ID;  Worthington,  OH;  accuracy  ± 
20  g)  before  and  after  each  experiment.  Water  from  premeasured 
bottles  was  available  to  drink  at  will  during  all  experiments,  and  a 
small  meal  (~500  kcal)  was  provided  during  the  8-h  experiments.  The 
weights  of  all  food  and  water  consumed  and  urine  voided  were 
measured  on  an  electronic  scale  (Ohaus  E1M210;  Nanikon,  Switzer¬ 
land;  accuracy  ±  1  g).  The  weight  of  any  fecal  mass  was  determined 
from  body  mass  changes  before  and  after  void.  Total  sweat  losses 


were  determined  from  mass  balance  using  the  Peters-Passmore  equa¬ 
tion  (3)  and  then  time  weighted  to  derive  rate.  The  measurements 
assume  that  1  ml  of  sweat  is  equal  to  a  mass  of  1  g.  To  summarize, 
sweat  loss  (kg)  =  change  in  body  mass  +  (solids  in  —  solids  out)  + 
(fluids  in  —  fluids  out)  —  (gases  in  —  gases  out),  where  gases  refer  to 
CO2-O2  exchange. 

These  data  sets  are  summarized  briefly  below.  Data  sets  l— IV  were 
used  to  develop  the  new  and  corrected  sweating  rate  prediction 
equations,  whereas  data  sets  V  and  V7  were  used  for  cross-validation. 

Data  set  I.  Details  can  be  found  in  Montain  et  al.  (15).  Data  were 
obtained  from  19  previously  heat-acclimatized  individuals  (13  men 
and  6  women)  who  completed  all  experiments  dressed  in  hot  weather 
battle  dress  uniforms  [BDU;  with  a  clo  (clothing  insulation  coeffi¬ 
cient)  value  of  1.08  and  a  im/clo  (evaporative  impedance  coefficient) 
value  of  0.49],  Volunteers  completed  12  randomized  exercise-heat 
stress  trials  in  which  they  walked  at  three  exercise  intensities  charac¬ 
terized  as  easy,  moderate,  and  heavy  work  intensities  set  at  250,  425, 
and  600  W  in  three  humid  environments  (Ta  =  28°C/PW  =  2.8  kPa, 
Ta  =  32°C/PW  =  3.57  kPa,  and  Ta  =  36°C/Pw  =  4.45  kPa,  where  Ta 
is  ambient  temperature  and  Pw  is  water  vapor  pressure).  In  the  other 
three  heat  stress  trials,  volunteers  walked  at  425  W  in  three  dry 
environments  (Ta  =  36°C/P„  =  1.47  kPa,  Ta  =  41°C/Pw  =  1.95  kPa, 
and  Ta  =  46°C/PW  =  2.53  kPa).  Dry  heat  stress  trials  were  completed 
following  a  humid  test  condition.  Appropriate  work-rest  cycles  for 
each  exercise  task  were  initially  determined  using  a  current  model 
(19)  for  a  predicted  2-h  total  exposure.  In  the  minority  of  heat/exercise 
trials  where  mean  skin  temperature  (TSk)  was  not  available,  TSk  was 
estimated  using  Saltin's  (20)  equation:  TSk  =  0.215Ta  +  26.6  (  ±  0.5 
SEE),  where  SEE  is  standard  error  estimate,  to  calculate  the  pertinent 
heat  balance  equation  parameters  essential  for  model  algorithms. 

Data  set  II.  Details  can  be  found  in  Cheuvront  et  al.  (3).  Thirty-nine 
healthy  volunteers  participated  in  this  study  specifically  designed  to 
obtain  data  for  the  development  of  an  improved  sweat  rate  prediction 
equation.  The  clothing  ensemble  was  the  U.S.  Army  woodland  BDU 
with  field  cap,  sleeves  down  [clo  =  1 .08,  im/clo  =  0.49  at  wind  speed 
( V)  =  1  m/s],  and  athletic  shoes.  Test  sessions  lasted  either  2  or  8  h. 
Twenty-one  volunteers  (16  men  and  5  women)  participated  in  the  2-h 
trials,  and  physical  characteristics  for  this  group  are  shown  in  Table  1. 
Eighteen  different  volunteers  (17  men  and  1  woman)  completed  the 
8-h  trials,  and  their  characteristics  are  shown  in  Table  2.  Tables  1  and 
2  also  describe  the  seven  different  levels  of  environmental  stress, 
work-rest  cycles,  metabolic  rate  at  each  exercise  intensity,  and  mea¬ 
sured  sweat  loss.  In  the  2-h  trials,  volunteers  were  not  heat  acclimated, 
whereas  in  the  8-h  trials,  volunteers  were  heat  acclimated. 

Data  set  III.  Details  can  be  found  in  Cheuvront  et  al.  (2).  Thirteen 
men  participated  in  this  study  specifically  designed  to  obtain  data  for 
the  development  of  an  improved  sweating  rate  prediction  equation. 
Each  subject  completed  three  trials  consisting  of  4  h  of  treadmill 
walking  (~500  W)  in  a  hot,  dry  environment  (Ta  =  35°C,  Pw  =  1.7 
kPa,  V  =  1  m/s).  The  U.S.  Army  BDU  was  worn  in  all  three  trials: 
alone  (trial  BDU:  clo  =  1.12,  im/clo  =  0.44  at  V  =  1  m/s),  combined 
with  interceptor  body  armor  (trial  IB  A:  clo  =  1.35,  im/clo  =  0.27  at 
V  —  1  m/s),  or  combined  with  IB  A  and  a  spacer  vest  (trial  SP:  clo  =  1.28, 
im/clo  =  0.32  at  V  =  1  m/s).  In  the  BDU  trial,  the  BDU  was  worn  with 
field  cap,  sleeves  down,  and  athletic  shoes.  The  IBA  vest  included  front 
and  rear  ballistic  protective  inserts  (throat  and  groin  protection  excluded). 
The  total  weight  of  the  vest  as  used  was  7.5  kg,  and  it  covered  —25%  of 
the  total  BSA.  In  trial  SP.  the  spacer  was  a  1-cm-thick  vest  of  proprietary 
knit  fabric  worn  between  the  IBA  and  uniform.  The  spacer  vest  is 
designed  to  produce  an  air  channel  that  theoretically  increases  the  poten¬ 
tial  for  ventilation  and  evaporative  cooling  of  the  torso. 

Data  set  IV.  Details  can  be  found  in  Chinevere  et  al.  (4).  One 
woman  and  five  men  did  continuous  treadmill  exercise  (—400  W)  for 
2  h  while  dressed  in  BDU  plus  IBA,  as  in  data  set  III  (IBA:  clo  = 
1.35,  im/clo  =  0.27  at  V  =  1  m/s).  Environmental  conditions  were  Ta  = 
307PW  =  2.1  kPa,  Ta  =  35°C/P„  =  4.27  kPa,  andTa  =  40°C/PW  =  1.47 
kPa.  All  other  procedures  and  details  were  as  outlined  in  Data  set  III. 
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Table  1.  Key  descriptive  and  physiological  data  for  2-h  experiments  in  data  set  II-B  used  in  developing  the  present  equation 

Trial 

Ta,  °C 

Pw,  kPa 

Work:Rest  Cycles,  #(min) 

n 

BSA,  m2 

Body  Weight,  kg 

TSk,  °C 

Vo2,  1/min 

SR,  1/h 

Al 

15 

0.85 

2X  (50:10) 

12M 

1.94±0.15 

77.6±8.8 

29.3±  1.0 

0.94±0.14 

0.1 35  ±0.079 

4F 

1.70±0.03 

64.3  ±4.2 

28.2±0.82 

0.79±0.17 

0.1 88  ±0.086 

A 

15 

0.85 

2X  (50:10) 

13M 

1.96±0.14 

79.7±9.3 

28.7±0.8 

1.49±0.17 

0.305±0.136 

2F 

1.68 

60.4 

28.7 

1.3 

0.319 

B 

15 

0.85 

2X  (50:10) 

13M 

1.96±0.14 

79.7±9.3 

29.1  ±1.0 

1 .99  ±0.2 1 

0.472  ±0.1 70 

2F 

1.68 

6.4 

27.5 

1.57 

0.424 

C 

20 

1.17 

2X  (50:10) 

14M 

1.96±0.14 

80.7 

31.1  ±0.8 

1 .08  ±0. 15 

0.220±0.089 

IF 

1.63 

56.9 

31.34 

0.81 

0.165 

D 

20 

1.17 

2X  (50:10) 

15M 

1.98±0.13 

80.9±8.9 

30.2±  1.08 

1.45±0.19 

0.410±0.177 

4F 

1.67±0.04 

62.6±5.7 

30.4  ±0.47 

1 . 14  ±0. 13 

0.346±0.084 

E 

20 

1.17 

2X  (50:10) 

15M 

1.97±0.14 

80.5  ±9.4 

30.3±  1.15 

2.03  ±0.26 

0.625  ±0.199 

2F 

1.69 

63.5 

29.6 

1.49 

0.498 

F 

25 

1.59 

2X  (50:10) 

11M 

1.96±0.11 

80.2±8.9 

3 1.8  ±0.55 

1 .03  ±0. 16 

0.321±0.102 

2F 

1.68 

60.35 

31.1 

0.84 

0.421 

G 

25 

1.59 

2X  (50:10) 

9M 

1 .95  ± .  12 

80.0±9.5 

3 1.5  ±0.92 

1.47  ±0.20 

0.479±0.171 

IF 

1.73 

63.8 

30.9 

1.26 

0.515 

H 

25 

1.59 

2X  (50:10) 

10M 

1 .95  ±0. 1 1 

80.5±8.9 

32.1  ±0.62 

1.99±0.27 

0.755±0.237 

IF 

1.63 

56.9 

31.6 

1.54 

0.558 

I 

30 

2.12 

2X  (50:10) 

10M 

1 .95  ±0. 1 1 

80.5±8.9 

33.5±0.75 

1.98  ±0.26 

0.935±0.317 

IF 

1.63 

56.9 

33.1 

1.45 

0.702 

Values  are  means  ±  SD  except  when  n  <  3;  n  =  no.  of  subjects  (M,  male;  F,  female)  completing  the  experiment.  Ta,  air  temperature;  Pw,  ambient  water  vapor 
pressure;  TSk,  mean  skin  temperature;  SR,  observed  sweating  rate.  Wind  speed  for  all  experiments  was  1  m/s.  Trial  letters  are  as  presented  in  Cheuvront  et  al.  (3). 


Data  set  V.  An  archival  raw  data  set  was  obtained  from  a  US- 
ARIEM  field  study  conducted  at  Ft.  Bliss,  TX,  to  determine  sweating 
rates  during  walking  exercise  in  mild  to  moderate  solar  load  condi¬ 
tions.  Details  can  be  found  in  Santee  et  al.  (21).  In  brief,  eight  males 
[average  body  mass  =  80.5  ±  15.2  kg  (SD);  BSA  =  1.97  ±  0.18  m2 
(SD)]  walked  at  2  miles/h  for  12  miles  (24  min  of  continuous  exercise 
with  a  6-min  break)  on  a  calibrated  track  for  6  h.  Subjects  carried  a 
22-kg  pack,  and  average  heat  production  was  maintained  at  ~350  W. 
Standard  chemical  protective  clothing  designated  as  mission-oriented 
protective  posture  (MOPP)  was  worn  during  the  walks.  The  clothing 
system  configurations  consisted  of  MOPP  0  (BDU:  clo  =  1.34  and 
im/clo  =  0.31  at  V  =  1  m/s),  MOPP  1  (clo  =  1.97  and  im/clo  =  0.17 
at  V  =  1  m/s),  and  MOPP  4  (clo  =  2.44  and  im/clo  =  0.12). 
Volunteers  walked  at  night  and  day  in  the  various  MOPP  clothing 
configurations.  Daytime  ambient  conditions  were  Ta  =  23°C/PW  = 
0.2  kPa,  with  a  mixed  solar  load  (diffuse  and  direct)  estimated  by  the 
effective  radiant  field  (ERF)  at  500  W/m2,  resulting  in  an  operative 
temperature  (T0)  of  49.5°C  (1,  5,  9,  27). 

Data  set  VI.  Raw  data  from  an  international  cooperative  data 
sharing  program  (8)  were  used  to  test  the  developed  equations.  Data 
from  13  males  and  9  females  (the  latter  in  the  follicular  phases  of  their 
menstrual  cycle)  were  compiled  to  compare  the  various  equations. 
Physical  characteristics  (±SD)  were  body  mass  =  82.7  ±  12.5  kg  and 
BSA  =  2.01  ±  0.16  m2  for  males  and  body  mass  =  60.4  ±  8.9  and 


BSA  =  1.66  ±  0.15  m2  females.  The  ambient  conditions  were  Ta  = 
40°C  and  Pw  =  1.47  kPa  at  V  =  0.4  m/s.  Subjects  exercised  for  2  h 
(4  exercise-rest  cycles  of  15:15  min)  or  until  their  rectal  temperature 
reached  a  peak  value  of  39.5°C  or  heart  rate  became  elevated  higher 
than  180  beats/min  for  3  min.  Average  metabolic  rate  was  ~365  and 
— 338  W  for  the  male  and  female  groups,  respectively.  Subjects  were 
dressed  in  the  Canada  nuclear,  biological,  chemical  (NBC)  protective 
clothing  system  previously  evaluated  using  USARIEM  manikin  pro¬ 
cedures  and  found  to  have  heat  transfer  characteristics  of  clo  =  1.88 
and  im/clo  =  0.18  at  V  —  1  m/s  (1). 

Heat  transfer  analyses.  Each  element  of  the  comprehensive  heat 
balance  equation  (5)  was  determined  from  the  raw  data  and  concate¬ 
nated  together  in  a  unified  spreadsheet  (Microsoft  Excel)  for  later 
analysis.  The  techniques  documented  by  Gagge  and  Gonzalez  (5)  and 
Breckenridge  (1)  were  applied  to  the  raw  spreadsheet  values  to 
determine  respective  heat  and  mass  transfer  coefficients  for  calculat¬ 
ing  dry  heat  exchange  ( R  +  C),  Emax,  and  other  variables  in  the  heat 
balance  equation.  Clothing  heat  and  evaporative  potential  parameters 
were  determined  on  a  regionally  heated,  articulated  manikin  at  various 
wind  speeds.  A  heat  balance  analysis  was  carried  out  for  each 
individual  response. 

In  general,  heat  balance  (in  W/m2)  of  the  human  body  surface  area 
can  be  expressed  by 


Table  2.  Key  descriptive  and  physiological  data  for  8-h  experiments  in  data  set  II-B  used  in  developing  the  present  equation 


Trial 

Ta,  °C 

Pw,  kPa 

Work:Rest  Cycles,  #(min) 

n 

BSA,  m 2 

Body  Weight,  kg 

TSk,  °C 

Vo2,  1/min 

SR,  1/h 

j 

40 

2.95 

6X  (60:20) 

12M 

2.03  ±0.14 

84. 7±  12.4 

35.9±0.61 

1.11  ±0.21 

0.667±0.114 

IF 

1.63 

55.3 

36.2 

0.775 

0.404 

K 

35 

1.69 

6X  (60:20) 

15M 

1 .97  ±0. 12 

80.2±9.8 

34.4±0.5 

1.39±0.21 

0.569±0.062 

IF 

1.62 

54.8 

34.12 

0.92 

0.405 

L 

35 

1.69 

6X  (60:10) 

15M 

1 .99  ±0. 15 

8 1 .7  ±  12.9 

33.9±0.5 

1 .05  ±0. 15 

0.452  ±0.058 

M 

27 

1.43 

6X  (60:20) 

14M 

1 .98  ±0. 15 

8 1 .7  ±  12.9 

32.3  ±0.6 

1.44  ±0.23 

0.406  ±0.096 

IF 

1.63 

55.7 

31.83 

0.96 

0.263 

N 

27 

1.43 

6X  (60:20) 

12M 

1 .98  ±0. 16 

81.1  ±  13.9 

32.7  ±0.64 

1.05±0.14 

0.269±0.05 

IF 

1.63 

55.6 

33.0 

0.770 

0.153 

O 

20 

1.17 

6X  (60:20) 

12M 

1.97±0.16 

80.9  ±14 

30.6  ±  1.15 

1.39±0.22 

0.229±0.087 

IF 

1.63 

55.1 

30.3 

0.847 

0.162 

Values  are  means  ±  SD  except  when  n  <  3;  n  =  no.  of  subjects  completing  the  experiment.  All  subjects  were  previously  heat  acclimated.  Wind  speed  for 
all  experiments  was  1  m/s.  Trial  letters  are  as  presented  in  Cheuvront  et  al.  (3). 
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PREDICTION  EQUATIONS  OF  SWEAT  LOSS 


S  =  net  metabolic  heat  flux  —  skin  insensible  heat  flux 

(2) 

—  skin  sensible  heat  loss, 

S  =  (M  ±  Wk)  -  E  ±  (R  +  C  +  K)  (3) 

where  S  is  the  time  rate  of  change  of  body  heat  (gain  or  loss).  If 
positive,  the  body  is  increasing  its  core  and  skin  temperature,  and 
these  can  be  estimated  using  M,  the  rate  of  metabolic  heat  production; 
Wk,  the  rate  of  accomplished  mechanical  work;  E,  the  rate  of  evap¬ 
orative  heat  loss  via  regulatory  sweating  from  eccrine  sweat  glands, 
diffusion  (Edif),  respiration  (£res),  and  metabolic  heat  loss  (mr);  C,  the 
rate  of  convective  heat  loss  from  the  total  body  surface  and  respira¬ 
tion;  R,  the  rate  of  radiant  heat  loss  (or  gain  from)  the  surrounding 
surfaces,  and  K,  the  rate  of  conductive  heat  flux  to  or  from  the 
environment. 

Radiation  exchange.  In  any  thermal  environment,  a  linear  radiation 
transfer  coefficient  (hr,  in  W-m_2-°C_1)  may  be  derived  (1,  5)  by 

hr  =  4ao-(Ar/AD)  •  fcl •  (5 .67  X  1(T8)  [(T„  +  Tsurf)/2  +  273.15]3  ( 4 ) 

where  a  is  the  skin  or  clothing  absorptance  for  the  radiation  exchange 
to  the  ambient,  a  is  the  Stefan-Boltzmann  constant  (5.67  X  10  ~s,  in 
W-m~2-K~4),  and  the  factor  A,/AD  is  the  ratio  of  the  effective 
radiating  area  of  the  human  body  to  the  total  body  surface  area  as 
measured  by  the  Dubois  surface  area  formula  (0.72  for  standing 
individuals).  The  interior  environmental  temperature  is  composed  of 
an  average  of  the  operative  temperature  (T0)  plus  all  the  surface 
temperatures  (Tsurt)  including  any  clothing  surface  temperature  (Tci); 
fci  represents  a  factor  that  increases  the  effective  Ar  of  the  body 
surface  by  some  15-20%  per  clo  (1,  5).  Mean  clothing  surface  is 
derived  from  Tsk  in  the  relationship  [T„  +  Fci(Tsk  —  T0)].  For  shorts 
and  a  T-shirt,  Fci  =  1  and  Tci  =  Tsk. 

Convective  exchange.  This  heat  exchange  factor  is  represented  by 
free  convection  and  forced  convection  via  increased  metabolic  activ¬ 
ity  or  increased  room  air  movement  artificially.  Two  equations  for 
estimating  the  convective  heat  transfer  coefficient  (hc,  in  W-m~2,°C) 
have  been  formulated  based  on  a  composite  of  free  and  forced 
convection  (1,  5,  9): 

hc  =  1.2  [(M  -  50)(PB/760)]°39  (5) 

where  M  is  metabolic  heat  production  (in  W/m2)  and  PB  is  the 
barometric  pressure  in  Torr  (1  kPa  is  equivalent  to  7.5  Torr  at  sea 
level).  Alternatively,  hc  for  fan-generated  forced  convection,  in  which 
ambient  air  movement  (V,  m/s)  is  the  main  factor  affecting  convective 
heat  exchange,  can  be  expressed  by  either 

hc  =  8.6  [V-PB/760]°53,  (6) 

when  persons  are  dressed  in  shorts  and  a  T-shirt,  or  by 

hc  =  12.7  [V-PB/760]°'5°,  ( 6 ') 

when  persons  are  clothed  (5). 

Insensible  heat  loss.  Procedures  applicable  to  clothing  heat  transfer 
were  used  to  calculate  environmental  heat  exchange  and  applied  to  the 
respective  data  for  each  individual  (1,  5,  6,  9).  These  methods 
consider  the  skin,  clothing,  and  environment  as  a  total  system  and  the 
constants  defining  insulation  and  water  vapor  transfer  as  functions  of 
effective  air  movement  (Veff).  The  term  Veff  is  the  sum  of  air  motion 
around  a  stationary  object  plus  the  speed  at  which  the  object  is 
moving.  IT(Veff)  is  the  total  resistance  to  heat  flow  by  radiation  and 
convection  (in  clo  units,  1  clo  is  equivalent  to  0.155  nr-K-W-1  or 
thermal  conductance  of  6.45  W-m_2-0C_1),  and  im(Veff)  is  the 
relative  total  resistance  to  evaporative  heat  transfer  (zero  to  one, 
dimensionless).  In  heat  balance  calculations,  im  is  not  used  alone  but 
as  a  latent  heat  transfer  coefficient  (im/IT),  separately  evaluated  using 
an  articulated,  moving,  sweating  manikin;  this  latter  quantity  is 
considered  as  a  key  dynamic  constant  incorporating  both  heat  and 
mass  transfer  via  “pumping”  through  cuffs,  vents,  and  walking; 


enthalpic  evaporation  (7,  10);  and  relative  permeation  from  skin  to 
each  subsequent  intrinsic  clothing  layer,  and  ultimately  to  the  ambient 
temperature  important  in  total  latent  heat  transfer  efficiency  of  mili¬ 
tary  clothing  present  in  this  study.  These  factors  were  not  considered 
in  development  of  the  legacy  Em ax  of  the  Shapiro  et  al.  (28)  study  and 
will  likely  improve  the  efficacy  of  the  OSE. 

Values  for  IT  and  (im/Tr)  as  a  function  of  Veff  were  calculated  for 
each  clothing  system  in  each  calculation  of  £jeq  in  the  heat  balance 
output  during  each  exercise  condition  or  use  of  body  armor  and 
transient  dynamic  effects.  These  clothing  parameters  were  ascertained 
from  the  following  power  curves  automatically  estimated  on  the 
sweating,  articulated  manikin  used  to  evaluate  clothing  ensembles  (1): 

lT  =  A-(Veff)B  (7) 

and 

(im/IT)  =  C-(Veff)D  (S) 

where  the  coefficients  A  and  C  are  the  values  for  IT  and  (im/IT)  when 
Veff  =1.0  m/s,  and  the  coefficients  B  and  D  are  slopes  of  plots  of  In 
(It)  and  In  (im/IT)  vs.  In  (V). 

The  intrinsic  thermal  insulation  value,  Iini(V),  was  obtained  by 
subtracting  the  value  of  the  insulation  of  the  air  boundary  layer,  Iaci, 
from  IT: 


1 

Iint  =  It - .  (9) 

(fac,M0.61  +  1.87^) 

where  faci  in  Eq.  9  is  the  increase  in  surface  area  due  to  clothing  that 
is  estimated  (1,  5)  using 

fad  =  1  +  (0.2- AD).  (10) 


The  algebraic  sum  of  the  total  (dry)  heat  loss  by  radiant  and  convec¬ 
tive  heat  exchange  ( R  +C),  in  watts  (1,  9),  therefore,  is 


Dry 


(R  +  C) 


6.45  •  ADfacl 

Ij  lit 


0-61(fsk  -  ir)  +  #-(T,k— TJ 
.'  +U0.6I  +  1.87  ^V) 

-Mnt 


(W 


Insensible  heat  loss  (£j  was  determined  by  the  rate  of  sweat 
secretion  (msw)  and  the  maximal  rate  of  evaporative  heat  loss  from  a 
fully  wetted  skin  surface  (£max).  £max  is  a  function  of  the  vapor 
pressure  gradient  between  the  fully  wetted  skin  surface  and  the  air 
(Ps.sk  —  Pw),  the  evaporative  heat  transfer  coefficient  (he),  and  im, 
Woodcock’s  dimensionless  factor  for  permeability  of  water  vapor 
through  clothing  (1,5,  6).  The  evaporative  heat  transfer  coefficient,  he, 
is  directly  related  to  the  convective  heat  transfer  coefficient,  he,  by  the 
Lewis  relationship  (LR;  2.2°C/Torr  or  16.5  K/kPa)  (5,  7,  10). 

When  evaporation  is  not  restricted  by  clothing  or  the  environment, 
then 


E  sk  =  ms„-\  (12) 

where  msw  is  in  grams  per  hour  and  A.  is  the  heat  of  vaporization  for 
sweat  at  35°C  (0.68  W-h-g-1),  as  determined  by  Wenger  (31).  The 
expression  for  skin  evaporative  heat  loss  (£Sk)  under  conditions  where 
evaporation  of  sweat  is  restricted,  particularly  during  exercise  with 
impermeable  clothing  (1,  5,  7,  10)  and  where  there  is  frank  dripping 
(Edrip)  or  wasted  sweat  due  to  skin  wettedness  (go  >  1.0)  (5,  7),  is 

£sk  =  (0.06  +  0.94oo)-AD£max  =  (LR-6.45)AD-(in/IT)-(Ps,sk-Pw), 

(12') 

when  £max  £  -msw-\,  where  AD  is  the  DuBois  surface  area  (m2)  (1, 
5)  and  Ps,sk  (in  Torr)  is  the  vapor  pressure  of  saturated  air  at  skin 
temperature.  Ps,sk  is  related  to  Tsk  by  the  Antoine  equation  (5,  7): 
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{  4030.183  \ 

Ps,sk  =  exP(l8.6686  -T-^235j  («) 

Respiratory  heat  loss  is  also  a  part  of  the  NSL  avenue.  Cres  and  Ercs 
are  directly  related  to  ventilation  rate  and  vary  as  a  function  of  aerobic 
exercise  intensity  (Mto t)  up  to  maximal  levels.  The  combined  equation 
for  estimating  respiratory  loss  by  convection  and  evaporation  (CTes  + 
£res,  g/min)  was  taken  from  Mitchell  et  al.  (18),  applicable  for  high 
levels  of  exercise  as 

Cres  +  Eres  =  0.019- Vo2(44  ~  PJ  (14) 

or 

(Cres  +  £res)  =  A d-Mm  [0.0014- (34  -  Ta)  +  0.0023  (44  -  Pj] 

(14') 

£res  is  also  modified  by  a  constant  (F)  for  high  levels  of  exercise  (20) 
so  that  if  V02  <  2.6  1/min,  F  =1,  and  if  V02  >  2.6,  F  =  1  + 
0.106(Vo2  -  2.6)2. 

Evaporative  heat  loss  corrected  for  metabolic  heat  losses  by  CO2 
and  O2  (mres,  g/min)  was  calculated  using  the  following  assumptions 
(20).  If  the  respiratory  exchange  ratio  (R)  is  equal  to  1, 

hires  =  V0,(R-pC02  —  pO,)  (15) 

where  pCCL,  the  density  of  CO2,  =  1.96  g/1  STPD;  p02,  the  density 
of  O2,  =  1.43  g/1  STPD,  and  mres  =  0.53(Vo2).  If  R  >  1, 

hires  =  Vo2(R-0.53)  U5') 

Statistical  procedures  and  quasi-Newton  analyses.  All  concate¬ 
nated  data  from  the  spreadsheet  were  analyzed  using  fuzzy  piecewise 
linear  and  nonlinear  regression  analyses  (29,  30)  with  various  statis¬ 
tical  modules  (Statistica,  version  7;  Tulsa,  OK)  to  establish  appropri¬ 
ate  change  points  in  sweat  loss  per  time  points  in  the  data  set,  coded 
by  trial,  sex,  and  individual  subject  number,  and  to  obtain  intercepts 
for  independent  parameters  derived  from  the  heat  balance  equation 
(£Veq,  £mai  )•  A  quasi-Newton  method  was  employed  to  derive 
regression  parameters.  In  this  method,  the  slope  of  a  function  at  a 
particular  locus  is  computed  as  the  first-order  derivative  of  the 
function  (at  that  locus).  The  "slope  of  the  slope”  is  the  second-order 
derivative,  which  documents  how  fast  the  slope  is  changing  at  the 
respective  point  and  in  which  direction.  The  quasi-Newton  method 
will,  at  each  step,  evaluate  the  function  at  different  points  to  estimate 
the  first-order  and  second-order  derivatives.  It  will  then  use  this 
information  to  follow  a  path  toward  the  minimum  of  the  loss  function 
(29).  The  fuzzy  piecewise  routine  is  more  robust  than  conventional 
methods,  is  not  sensitive  to  outliers,  and  tracks  transient  responses  or 
irregular  data  especially  found  in  this  data  set  due  to  heavy-intensity 
exercise  and  disparate  sweating  patterns.  The  technique,  as  con¬ 
structed  for  this  study,  is  also  suitable  for  long-term  time  series 
predictions  of  specific  variables  (29). 

OSE-predicted  sweating  rates  were  initially  compared  against  ob¬ 
served  data  to  obtain  residual  analyses  to  ascertain  how  much  OSE 
deviated  from  the  observed  data.  Next,  the  new  fuzzy  piecewise  (PW) 
equation  (29)  was  compared  with  OSE  and  the  observed  raw  data 
secured  for  each  separate  trial.  Corrections  to  OSE  (designated  as 
OSE,c)  were  derived  by  independent  piecewise  regression  analyses 
incorporating  an  iterative  approach  to  obtain  the  most  optimum 
equation  (exponential,  log  fit)  and  test  the  significance  of  the  derived 
regression  coefficients  (Wald  statistic)  that  fit  the  database  (17).  The 
conditions  were  that  the  standard  error  estimate  (SEE)  of  a  new 
equation  (PW)  should  not  deviate  from  the  observed  data  or  the 
independently  determined  fuzzy  piecewise  equation  by  more  than  ±125 
g-m~2-h~‘  (roughly  less  than  ±0.24  1/h  for  a  person  with  1.9  m2 
BSA).  This  is  a  more  liberal  criterion  than  the  one  used  previously  by 
Cheuvront  et  al.  (3),  which  assigned  an  a  priori  “zone  of  indifference” 
of  ±0.125  1/h  (i.e.,  ~65.8  g-m_2-h_1  for  a  1.9  m2  BSA)  resulting  in 


band  differentiation  of  some  ±2  g/min  over  an  extended  period.  The 
higher  level  relaxes  the  prediction  of  Ereq  that  can  be  calculated  from 
the  heat  balance  equation.  This  decision  is  based  on  the  fact  that  the 
calculation  of  evaporative  heat  exchange  parameters  in  clothed  indi¬ 
viduals  are  not  exact,  and  indeed,  the  latent  heat  of  evaporation  is 
variable,  particularly  during  cold  ambient  conditions  (10),  and  a 
model’s  accuracy  predicting  sweating  rate  is  not  always  tied  in  with 
evaporative  cooling  efficiency  (7,  9,  10),  particularly  when  w  is  100% 
for  an  extended  period. 

Data  are  expressed  as  means  ±  SD,  means  ±  SE,  or  means  ±  95% 
confidence  interval  (Cl).  The  differences  in  observed  sweating  rate, 
heat  production,  and  the  output  from  the  various  prediction  equations 
were  analyzed  using  a  factorial  ANOVA  design  to  include  main 
effects  and  interactions  for  categorical  predictors  (sex,  all  trials).  Both 
univariate  (using  a  given  single  continuous  dependent  variable)  and 
multivariate  designs  (multiple  continuous  dependent  variables)  were 
analyzed.  If  a  significant  F  value  was  found  for  a  given  dependent 
variable,  the  more  conservative  Bonferroni  adjustment  procedure  (17) 
was  used  as  a  post  hoc  approach  to  locate  critical  differences  at  P  < 
0.01  and  considered  statistically  significant;  correlations  with  a  P 
value  >0.01  (including  those  with  P  values  between  0.01  and  0.05) 
were  considered  nonsignificant. 

We  determined  the  prediction  accuracy  of  a  given  equation,  com¬ 
pared  with  observed  sweating  rates  from  composite  trial  data,  using 
analyses  of  differences  of  mean  error  rates  as  described  by  Lim  et  al. 
(13).  For  each  trial  data  set,  the  algorithm  with  the  lowest  error  rate 
compared  with  the  observed  sweating  rate  is  assigned  rank  one,  the 
second  lowest  rank  two,  and  so  on,  with  average  ranks  assigned  in 
case  of  ties  as  formulated  by  Lim  et  al.  (13).  The  Friedman  ANOVA 
nonparametric  test  (17)  was  then  used  to  test  differences  in  mean  error 
rate  ranks  for  each  algorithm. 

Finally,  an  independent  cross-validation  analysis  of  the  fuzzy 
piecewise  equation  was  executed  against  two  independent  archival 
data  sets:  a  field  study  in  which  a  group  of  men  walked  in  NBC 
clothing  (at  MOPP  levels  0,  1,  and  4)  (19,  21)  and  a  laboratory  study 
in  which  men  and  women  walked  in  NBC  protective  clothing  con¬ 
ducted  at  DRDC  (8);  experimental  details  are  presented  above. 

RESULTS 

Table  3  provides  the  calculated  mean  heat  production  ( M , 
W/m2)  and  observed  (measured)  sweating  rates  (g-m_2-h_1) 
±SE  for  data  sets  I-IV  (methods)  separated  by  sex.  After 
ANOVA  tests,  a  Bonferroni  post  hoc  (17)  analysis  was  per¬ 
formed.  The  only  significant  difference  in  measured  sweating 
rates  between  men  and  women  occurred  with  the  moderate- 
intensity  exercise  trials  in  data  set  I  (P  <  0.01). 

Figure  1  plots  the  individual  sweating  rates  (Fig.  1A)  for 
data  sets  I-IV  and  the  calculated  residuals  (Fig.  IB)  used  to 
evaluate  OSE  and  develop  a  new  prediction  equation.  Notice¬ 
able  are  the  high  points  and  depressions  in  sweating  rates  (Fig. 
1A)  due  to  variation  in  exercise  intensity  and  environmental 
conditions  (particularly  cooler  environments)  among  data  sets 
I-IV.  The  residual  values  in  Fig.  IB  (comparing  measured 
minus  predicted  sweating  rates  for  each  data  point)  demon¬ 
strate  that  some  predicted  values  were  overestimated  by  100% 
or  more,  especially  during  higher  exercise  intensities  (data  set 
I),  and  underestimated  by  80%  during  cooler  trials  ( data  set  II). 
However,  for  many  of  the  data  sets,  OSE  was  within  ±20%  of 
observed  data  with  residuals  distributed  symmetrically  around 
the  zero  line,  particularly  during  low  exercise  intensities  and 
mild  heat  stress  conditions. 

The  data  set  was  next  examined  to  develop  optimum  regres¬ 
sion  parameters  that  would  satisfy  all  data  sets  sufficiently 
(within  a  SEE  of  ±125  g  •  m  2  •  h  1  criterion  and  coefficient  of 
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Table  3.  Heat  production  and  observed 

PREDICTION  EQUATIONS  OF  SWEAT  LOSS 

sweating  rate 

Description 

Data  Set 

Sex 

M,  W/m2 

Mean  ±SE 

OSR,  g-m“ 

Mean 

"2*h~> 

±SE 

n 

Easy  work 

I-L 

M 

163.04 

8.87 

295.21 

18.35 

32 

I-L 

F 

148.60 

12.18 

208.31 

25.17 

17 

Moderate  work 

I-M 

M 

224.00 

5.84 

407.13 

12.06* 

74 

I-M 

F 

252.24 

8.37 

332.59 

17.29* 

36 

Heavy  work 

I-H 

M 

329.18 

8.25 

443.96 

17.06 

37 

I-H 

F 

357.61 

11.83 

356.88 

24.46 

18 

BDU 

H-B 

M 

252.42 

3.53 

224.72 

7.30 

202 

n-B 

F 

231.82 

10.04 

210.19 

20.75 

25 

Moderate  work,  BDU 

III-B 

M 

284.66 

13.92 

353.21 

28.78 

13 

III-B 

F 

NA 

NA 

NA 

NA 

1 

Moderate  work,  BDU+  body  armor 

III-I 

M 

287.79 

13.92 

407.64 

28.78 

13 

Moderate  work,  BDU+  body  armor 

III-I 

F 

NA 

NA 

NA 

NA 

1 

Moderate  work,  BDU+  body  armor  with  vest 

III-S 

M 

289.52 

13.92 

420.47 

28.78 

13 

Moderate  work,  BDU+  body  armor  with  vest 

III-S 

F 

NA 

NA 

NA 

NA 

1 

Easy  work,  BDU+  body  armor  with  vest 

IV 

M 

213.07 

12.96 

362.12 

26.79 

15 

Easy  work,  BDU+  body  armor  with  vest 

IV 

F 

193.80 

28.98 

483.49 

59.91 

3 

Values  are  least-squares  means  between  males  and  females  for  combined  data,  followed  by  SE  values;  n  =  no.  of  observations.  M,  heat  production;  OSR, 
observed  sweating  rate;  BDU,  battle  dress  uniform;  NA,  not  estimated  due  to  small  sample  in  trial.  *P  <  0.01  (Bonferron;  post  hoc  test  within  trial  cells). 


determination  R 2  >  0.8)  as  shown  in  Fig.  2,  A-C.  Figure  2 A 
shows  that  OSE  accounted  for  74%  of  the  variance  but  had  a 
high  SEE  (181  g-m_2-h_1).  Since  the  data  were  collected  at 
different  times  and  contained  disparate  conditions,  a  nonlinear 
fuzzy  piecewise  regression  employing  a  quasi-Newton  solution 
provided  the  best  resolution  (29,  30).  Fuzzy  piecewise  regres¬ 
sion  using  the  Eieq  and  £max  values  as  input  parameters  pro¬ 
duced  the  following  equation  (PW): 

Sweating  rate  (g-m_2,h_1)  =  147  +  1 .527-Ereq  —  0.87 -£max 

(76) 

Figure  2B  provides  the  comparison  of  output  from  PW  with 
respect  to  observed  sweating  rates  for  data  sets  I— IV.  PW 
accounted  for  78%  of  the  variance  for  the  same  data  with  a 
much  smaller  SEE  (72  g-m_2,h_1).  For  this  analysis,  £j-eq  and 
Emax  were  determined  from  the  individual  heat  balance  equa¬ 
tion  by  using  new  heat  and  mass  transfer  coefficients  as 
explained  in  methods  and  were  transformed  to  sweating  rate  by 
division  of  W-nW2/(0.68  W-h-g^1)  (31). 

Since  it  was  determined  that  the  output  from  OSE  would 
probably  not  completely  unify  the  data  set  and  predict  ade¬ 
quately  because  of  the  imprecise  clothing  coefficients  deter¬ 
mined  from  a  static  manikin  and  the  low  wind  speeds  (10) 
applied  to  the  heat  balance  analyses  in  OSE  (28),  a  similar 
iterative  approach  was  run  to  correct  the  OSE  by  using  indi¬ 
vidual  analysis  of  £j-eq  and  7tmax  using  modern  analytic  ap¬ 
proaches  to  the  clothing  coefficients  in  the  data  sets  (methods). 
An  exponential  correction  to  Eq.  1  was  successfully  obtained 
after  various  statistical  algorithm  solutions  were  attempted. 
This  solution  produced  the  following  correction  applied  to  the 
OSE  equation  to  predict  sweating  rates: 

OSEc(g-m_2-h^1)  =  147-exp(0  0012'OSE)  (77) 

where  OSE  in  Eq.  77  is  the  uncorrected  output  from  the 
original  Eq.  1  (28). 

£req  and  £max  signify  the  new  values  obtained  with  the  new 
clothing  coefficients.  Figure  2 C  plots  the  comparison  of  output 
from  OSE,c  with  respect  to  observed  sweating  rates  for  all  data 
sets.  OSE,c  accounted  for  —64%  of  the  variance  but  with  a 


small  SEE  (±98.3  g-m^2^"1).  Noticeably,  the  regression 
coefficient  (0.603)  is  significantly  smaller  (P  <  0.01,  Wald 
statistic)  than  the  OSE  output  equation. 

Figure  3  provides  a  histogram  comparing  mean  outputs  for 
OSE,  OSEc,  and  PW  equations  against  measured  sweating 
rates  for  data  sets  I-TV.  The  OSE  demonstrated  the  widest 
variability  compared  with  the  measured  data  and  was  markedly 
elevated  for  the  moderate-  and  high-intensity  exercise  experi¬ 
ments.  However,  OSEc  and  PW  prediction  equations  tracked 
measured  sweating  rates  remarkably  well. 

Cross-validation  analyses.  To  evaluate  the  prediction  equa¬ 
tions  against  averaged  responses  from  independent  data  sets, 
we  made  comparisons  with  two  additional  data  sets:  one  an 
outdoor  field  study  ( data  set  V)  and  another  conducted  in  an 
environmental  chamber  ( data  set  VI).  In  both  studies,  volun¬ 
teers  wore  a  variety  of  clothing  configurations,  including  NBC 
protective  clothing  (see  methods).  Figure  4  plots  the  mean 
observed  sweating  rates  and  mean  predicted  values  during  a 
field  study  test  ( data  set  V)  (21).  In  general,  PW  and  OSEc 
equations  were  not  significantly  different  from  the  experimen¬ 
tally  observed  sweating  rates.  The  OSE,  however,  consistently 
predicted  responses  too  high  within  and  between  the  various 
trials  except  for  the  day  with  solar  load  and  night  trials  when 
subjects  wore  the  MOPP  level  1  clothing  configuration.  Figure  5 
plots  the  mean  observed  sweating  rates  and  mean  predicted 
values  during  the  environmental  chamber  study  ( data  set 
VI).  The  OSE  predicted  sweating  rates  higher  than  measured 
data  consistently  within  the  two  male  and  female  groups  of 
subjects  (P  <  0.01). 

DISCUSSION 

This  study  compared  the  accuracy  of  OSE  to  a  new  equation 
(PW)  for  predicting  sweating  rate  during  extended  exercise  (2 
to  8  h)  with  higher  metabolic  rates  and  contemporary  clothing 
systems  including  body  armor  and  developed  a  correction  to 
OSE  (OSEc)  that  is  applicable  and  easily  migrated  into  various 
existing  rational  and  operational  thermal  models.  The  hypoth¬ 
eses,  based  on  current  research  and  literature  results  (3,  15), 
were  that  OSE  overestimates  sweating  rates  and  that  improved 
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Fig.  1.  A:  individual  sweating  rate  measurements  from  data  sets  I-IV. 
B:  residual  plot  of  original  Shapiro  equation  output  (OSE)  (28)  vs.  observed 
data  for  each  individual  ( n  =  504). 


prediction  equations  could  be  developed.  The  newly  developed 
OSEc  and  PW  equations  were  better  predictors  of  sweating 
rates  (58  and  65%  more  accurate,  P  <  0.01)  and  produced 
minimal  error  (SEE  <  100  g-m_2-h  ')  for  conditions  both 
within  and  outside  the  original  OSE  domain  of  validity,  which 
include  cooler  environments,  higher  metabolic  rates,  longer 
work  durations,  and  modern  protective  clothing  and  equipment 
ensembles. 

The  rationale  to  develop  a  predictive  equation  for  sweating 
rate  (and,  thereby,  water  needs)  was  a  unique  concept  at  the 
time  that  OSE  was  developed  and  one  that  had  not  been 
expanded  on  until  this  effort.  The  OSE  predicts  sweat  losses, 
and  thus  water  needs,  over  wide  thermal  environments  by 
knowledge  of  only  two  key  variables,  Ereq  and  £max,  which 


Fig.  2.  Lineal'  plots  of  prediction  equation  output  vs.  observed  sweating  rates. 
A:  original  Shapiro  equation  (OSE).  B:  piecewise  equation  (PW).  C:  corrected 
OSE  (OSEc).  The  correlation  between  measured  sweating  rates  and  predicted 
sweating  rates  from  each  equation  is  indicated  in  each  panel.  Fuzzy  PW  is  a 
nonlinear  transient  analysis;  linear  correlations  are  shown  only  for  relative 
comparison  of  accuracy  between  the  various  output  algorithms.  SEE,  standard 
eiTor  of  the  estimate. 
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Fig.  3.  Sweating  rate  output  (means  ±  SD)  from  the  various  equations 
compared  with  measured  (observed,  Obs  SW)  data  plotted  for  data  sets  /-TV. 
Data  set  designations  are  the  same  as  in  Table  3.  Prediction  accuracy  and 
significance  was  determined  by  analyses  of  mean  error  rate  of  each  equation 
vs.  observed  data  (see  Table  4  and  text). 


directly  or  indirectly  integrate  the  effects  of  the  internal  factors 
(metabolic  rate,  skin  and  core  temperature)  and  external  factors 
(clothing,  operative  temperature,  wind,  and  humidity).  OSE 
has  been  shown  to  often  overestimate  sweating  rates  for  a 
variety  of  conditions  within  and  outside  the  original  equation’s 
domain  of  validity  (3). 

The  original  equation  was  not  meant  to  predict  responses  at 
high  metabolic  rates.  Elevated  heat  production  is  accompanied 
by  higher  levels  of  core  and  skin  temperatures  (20)  and 
generally  higher  sweating  rates  that  can  saturate  the  microcli- 
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Fig.  5.  Comparison  of  output  from  the  various  prediction  equations  and 
measured  (observed)  sweating  rate  plotted  for  each  group  (9  females  and  13 
males)  in  an  experimental  cross  validation  study  (study  VI).  Values  are 
means  ±  SD.  NS  differences  between  genders;*  =  P  <  0.01  within  group. 


mate  (and  invariably  Zsdrip),  both  of  which  become  limiting 
factors  in  a  person’s  ability  to  achieve  steady-state  heat  bal¬ 
ance.  Although  core  temperatures  are  not  shown  in  this  study, 
this  variable  is  indirectly  related  to  the  analysis  of  the  heat 


Fig.  4.  Comparison  of  output  from  each  of  the  prediction 
equations  and  measured  (observed)  sweating  rate  plotted  for 
each  group  walking  in  the  outdoor  cross-validation  study 
{data  set  V)  clothed  in  various  NBC  clothing.  MOPP,  mis¬ 
sion-oriented  protective  posture.  Values  are  means  ±  SD. 
*P  <  0.01  within  each  specific  trial,  a  =  P  <  0.01  between 
night  and  day  trials  in  MOPP  O. 
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balance  and  reflected  directly  by  any  changes  in  S.  The  latter  is 
affected  by  core  and  skin  temperature  drive  so  that  .S’/hFcl  = 
(dTb/dr)[(0.97mb/AD)]/hFcl,  where  Tb  is  integrated  mean 
body  temperature,  mb  is  body  mass  (kg),  AD  is  body  surface 
area  (m2),  0.97  is  body  specific  heat  (W-h)/(kg-°C),  and  hFcl 
is  the  combined  heat  transfer  coefficient  times  Burton’s  cloth¬ 
ing  factor,  Fcl  (5,  26). 

Also,  steady-state  skin  temperature  levels  add  errors  in 
determining  dry  heat  flux  (R  +  C  +  K)  through  clothing  and 
insensible  heat  loss,  and  variable  latent  heat  of  evaporation 
(10)  also  can  introduce  errors  in  £max  calculations  leading  to 
wide  discrepancies  in  prediction  accuracy  of  sweating  rates 
using  OSE. 

In  this  study,  OSE  consistently  overestimated  measured 
sweating  rate  (Fig.  2,  A-C),  particularly  at  the  higher  measured 
sweating  rates.  A  correction  to  OSE  using  the  present  extensive 
database  (101  subjects,  longer  work  durations  from  4  to  8  h, 
and  a  variety  of  clothing  systems)  extended  accurate  prediction 
of  actual  sweating  rates  from  —200  to  —500  g  •  rn  2  •  h  “ 1 .  The 
most  significant  finding  of  the  current  study  was  the  develop¬ 
ment  of  a  new  and  improved  equation  PW:  msw  (g-m  “2-h  ')  = 
147  +1.527  -  0.87 -£max. 

The  improved  equation  PW  is  based  on  fuzzy  piecewise 
regression  analysis  that  incorporates  the  combined  effects  of  a 
broader  range  of  metabolic  rates,  wearing  of  body  armor 
(commonly  employed  by  military  and  law  enforcement),  and 
other  aspects  of  modern  military  clothing  systems.  This  equa¬ 
tion  is  essentially  applicable  for  both  men  and  women  working 
for  time  periods  up  to  8  h.  Most  importantly,  this  equation  still 
incorporates  attributes  of  OSE  documenting  two  essential 
physiological  mechanisms  (5,  26)  necessary  in  thermoregula¬ 
tory  function,  one  coupled  with  eccrine  sweat  gland  function, 
itself  associated  with  internal  and  skin  body  temperatures 
necessary  in  the  regulation  of  heat  balance  (£req),  and  the  other 
needed  whenever  heat  transfer  mechanisms  must  be  character¬ 
ized  in  terms  of  clothing  worn  and  environmental  impact 
(£max)-  In  the  latter  cases,  the  factor  £j-eq/msw  is  a  direct 
correlate  of  M  (metabolic  heat  production),  since  Ereq  is  based 
on  the  solution  of  the  comprehensive  heat  balance  equation 
comprising  S  =  (M  ±  Wk)  — E±(R  +  C  +  K)  (10). 

The  relative  accuracy  of  PW  and  OSEc  has  only  been 
implied  so  far  in  this  discussion.  The  concern  is  that  any 
equation’s  output  should  accurately  predict  sweating  rates 
within  a  given  boundary  level  because,  over  an  8-h  period,  it 
has  an  impact  on  water  requirements  needed  by  an  individual. 
Cheuvront  et  al.  (3)  estimated  that  this  should  be  ±0.125  1/h 
(some  65.8  g-m~2-h~'  for  a  person  1.9  m2  BSA).  Our  results 
in  this  study  showed  via  the  SEE  (Fig.  2,  A-C)  that  the  various 
equations’  domains  of  accuracy  were  ±12 A  g  •  m  2  •  h  1  for 
PW,  ±98.3  g-m“2-h_1  for  OSEc,  and  ±181.3  g-m^-h'1 
using  the  original  OSE.  These  translate  for  the  1.9  m2  BSA 
person  to  ±0.137  1/h  (new),  0.186  1/h  (corrected),  and  0.340  1/h 
(original)  for  the  whole  database  considered  in  this  study.  It  is 
doubtful  that  the  ±0.137  1/h  found  using  the  new  equation  is 
physiologically  meaningfully  compared  with  the  ±0.125  1/h 
criterion,  which  would  result  in  —1,096  ml  over  8  h  com¬ 
pared  with  the  1,000-ml  criterion  recommended  by  Cheu¬ 
vront  et  al.  (3). 

A  sweating  rate  predictive  model  that  is  operational  for  wide 
uses  such  as  determining  water  requirements  over  various 
environments,  activities,  and  time  domains  is  only  as  good  as 


the  accuracy  it  possesses  compared  with  the  true  (or  experi¬ 
mentally  observed)  outcome.  Within  the  scope  of  this  study’s 
intent,  any  equation  developed  must  exhibit  a  tight  association 
between  the  model's  outcome  predictions  and  the  group  out¬ 
comes  of  the  test  population.  The  equation  serves  as  a  surro¬ 
gate  of  actual  observed  sweat  responses  based  on  one  or  more 
independent  factors  (Evcq,  £max)  driving  the  outcome  (mea¬ 
sured  msw).  The  analysis  of  mean  error  rates  for  all  combined 
trials  resulted  in  a  rejection  of  the  null  hypothesis  that  all  three 
equations  are  equally  accurate. 

Table  4  shows  that  PW  was  almost  65%  more  accurate  in 
predicting  sweating  rate  compared  with  the  OSE.  However, 
using  OSEc  would  still  be  58%  more  accurate  than  using  OSE. 
Therefore,  using  either  PW  or  OSEc  could  be  considered 
mathematically  equivalent  over  the  range  of  environments, 
time  domains,  or  metabolic  activities  studied.  Furthermore, 
integration  with  the  OSE  that  resides  in  many  rational  and 
operational  thermal  models  is  now  possible  by  applying  OSEc: 
msw  (g-m~2-h_1)  =  147-exp  (0.0012-OSE),  where  OSE  is 
output  from  OSE  (28)  that  leaves  the  exponential  regression 
coefficients  [27.9 -Ereq- (£max)-0'455]  intact  but  modifies  the 
Ereq  and  £max  to  conform  with  more  robust  and  modern  heat 
and  evaporative  transfer  coefficients  and  dynamic  manikin 
thermal  resistance  values. 

The  limits  of  the  corrected  predictive  equations  comprise  the 
original  £req  and  £max  limits  (50  <  Ereq  <  360  W/m2  and  20  < 
£max  <  525  W/m2,  respectively),  but  these  limits  are  extended 
so  that  the  equation  is  applicable  for  higher  metabolic  intensi¬ 
ties  ( M  =  400  W/m2,  —700-800  W),  lower  ambient  condi¬ 
tions  (Ta  =  15°C),  and  longer  time  periods  (up  to  8  h). 

PW  and  OSEc  were  derived  principally  from  environmental 
chamber  experiments.  The  limited  results  from  the  field  test 
reported  in  Fig.  4  ( data  set  V)  indicate  that  PW  and  OSEc 
appear  valid  in  predicting  sweating  rate  during  mild  effective 
radiant  loads.  Although  solar  load  in  terms  of  ERF  did  include 
daytime  tests  as  high  as  500  W/m2,  a  more  extensive  database 
incorporating  wider  thermal  conditions  with  more  intense  and 
variable  solar  loads  (27)  and  further  testing  or  modification  of 
the  equations  studied  is  needed  in  the  future. 

Conclusions.  The  present  study  developed  both  a  new  (PW) 
and  a  corrected  (OSEc)  sweating  rate  (and  fluid  needs)  predic¬ 
tion  equation  that  provides  a  more  accurate  estimate  of  mea¬ 
sured  sweating  rates,  and  therefore  also  water  needs,  for  a 
variety  of  military,  occupational,  and  sports  medicine  scenar¬ 
ios.  PW  and  OSEc  have  wide  application  to  public  health, 
military,  occupational  and  sports  medicine  communities.  Fu¬ 
ture  studies  should  develop  sweating  rate  databases  for  other 


Table  4.  Summary  of  mean  error  rate  analyses  for  each 
equation  compared  with  data  over  all  trials 


Variable 

Mean  ±  SD 

OSE  rank 

9.75  ±  2.3 1  * 

OSEc  rank 

5.37±3.50f 

PW  rank 

4.12±2.16 

OSE,  original  shapiro  equation;  OSEc,  corrected  OSE;  PW,  piecewise 
equation.  Rank  analysis  was  performed  using  Friedman  ANOVA.  *P  <  0.01, 
PW  vs.  OSE  (64.5%  difference).  fP  <  0.01,  OSEc  vs.  OSE  (57.9%  differ¬ 
ence).  There  is  no  significant  difference  comparing  PW  vs.  OSEc  (15.6% 
difference). 
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conditions  to  test  the  extension  validity  of  the  new  equation  for 
more  intense  and  variable  solar  load  conditions. 
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